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YIELD  OF  PHOSPHORUS  IN 
EASTERN  LAKE  ONTARIO  REGION 

Dr»  Syed  S.  Ahmed  Rizvi  1 
and  S.  E.  Salbach  2 


ABSTRACT 
The  yields  of  phosphorus  were  estimated  from  four  selected 
Eastern  Lake  Ontario  region  watersheds.    These  streams  drain  mainly 
agriculture  and  pasture  lands  to  Lake  Ontario.    The  variation  of  con- 
centrations, flows  and  yields  for  the  four  watersheds  were  calculated 
and  plotted  for  comparison  on  a  monthly  basis  over  the  four  year  study 
period.    Mathematical  equations  for  forecasting  of  concentration  levels 
were  developed  and  tested  for  statistical  significance.    The  mean 
monthly  yields  were  also  plotted  to  assess  on  a  quantitative  basis 
during  the  calendar  year,  the  variation  of  the  discharge  of  phosphorus 
to  Lake  Ontario.    About  40  to  5  0  per  cent  of  the  annual  yields  were 
found  to  be  released  during  the  months  of  February,  March  and  April. 
A  step-wise  regression  analysis  was  used  for  developing  the  relation- 
ship between  concentration,  yield  and  land  use  for  this  region.    Three  year 
annual  yields  of  phosphorus  from  the  four  watersheds  ranged  between 
237  to  810  lbs/mi^-yr.    The  methodology  presented  could  be  used  on 
other  water  quality  data. 


As  st.  Prof,  of  Eng.,   S.  G.  W.  Univ.,   Montreal  and  Consult.  Eng .  , 
WQSB,  OWRC,  Toronto,  Ontario.     September,   1968. 

2     P.  Eng.,  WQSB,  OWRC,  Toronto,  Ontario. 


INTRODUCTION 

In  recent  years,  the  Great  Lakes  have  received  widespread 
attention  from  public  and  private  groups  who  were  concerned  about 
the  increasing  evidence  that  through  direct  waste  discharges  and 
tributary  drainage,  the  Great  Lakes  are  becoming  excessively  en- 
riched with  nutrients  -  nitrogen  and  phosphorus.    This  report 
presents  the  findings  of  an  analysis  of  nutrient  data  for  four  streams 
tributary  to  Lake  Ontario  in  the  north-eastern  part  of  the  Province  of 
Ontario.    The  objectives  of  the  analysis  were  to  determine  the  yield 
of  total  phosphorus  from  forest  and  pasture  land  and  to  develop  a 
forecasting  equation  based  on  streamflows  and/or  other  independent 
variables.    The  nutrient  data  used  for  this  analysis  was  obtained  by 
the  Ontario  Water  Resources  Commission  as  part  of  a  regular  water 
quality  monitoring  program  of  lakes  and  streams  in  the  Province  of 
Ontario. 

The  region  under  study  is  drained  by  four  streams,  the  Moira 
River,  Salmon  River,  Napanee  River  and  Wilton  Creek  and  for  the  pur- 
pose of  this  report  is  referred  to  as  Eastern  Lake  Ontario  Region  (Figure 
1).    The  selected  area  has  essentially  similar  socio-economic,  physical, 
and  climatologicaJ    features. 


PHYSIOGRAPHY  AND  LAND  USE 

The  general  physiography  is  outlined  in  Figure  1.    The  region 
has  four  primary  soil  divisions  as  follows: 

In  the  northern  part  of  the  region  the  Canadian  Shield  pre- 
dominates.   It  is  composed  of  precambrian  granite  outcrops  but  glacial 
scouring  has  caused  numerous  small  deep  bank  lakes.     Shelf  depressions 
are  usually  filled  with  marsh  or  scrub,  and  vegetation  is  sparse.    There 
is  little  or  no  agriculture  in  the  area. 

Bordering  the  Canadian  Shield  to  the  south  is  the  till  moraine, 
an  area  of  rough  land.    The  underlying  bedrocks  of  the  till  moraine  are 
sedimentary  limestones  of  the  Black  River  group.    The  moraines  of  the 
area  are  characterized  by  angular  fragments  and  blocks  of  limestone. 
The  surface  is  extremely  rough  and  mainly  used  for  raising  cattle  and 
a  few  flocks  of  sheep.    Crossing  this  moraine  belt  are  the  Moira  tribu- 
taries which  cut  up  to  100  ft.  into  the  bedrock.    A  number  of  these 
valleys  are  blocked  by  glacial  drift;    thus  isolating  long  narrow  lakes 
and  swamps. 

The  Napanee  Plains  (clay  and  limestone  plains)  are  flat  to  un- 
dulating south-east  of  the  till  moraine  area.    While  the  soil  Is  only  a 
few  inches  deep  over  most  of  the  region  some  deeper  glacial  till  occurs 
in  the  stream  valleys  and  toward  the  north.    Alluvial  deposits  exist  in 


the  valleys.    There  are  also  a  few  scattered  drumlins  of  the  long  thin 
type.    The  greatest  elements  of  relief  are  provided  by  the  Salmon  and 
Napanee  rivers.    They  cut  the  rockplain  to  a  depth  of  50  to  100  feet. 
In  this  area  hay  and  oats  are  the  most  prevalent  crops  and  cattle  and 
poultry  are  in  notable  number. 

To  the  west,  north  of  Belleville  are  the  till  and  sand  plains. 
This  area  is  part  of  the  Iroquois  plain  which  extends  along  the  north 
shore  of  Lake  Ontario  into  the  Niagara  fruit  belt.    The  region  is  identi- 
fied by  cliffs,  boulder  pavements  and  shallow  soils.    Even  though 
rurally  populated  the  till  and  sand  plains  have  very  little  agriculture. 

Drumlins  exist  scattered  in  the  eastern  part  of  the  region, 
west  of  the  Moira  River.    This  is  composed  of  Grand  ridge  and  poor 
soils  (clays).    Along  the  northern  edge,  the  drumlins  are  small, 
scattered  and  imperfectly  formed.    Along  the  southern  edge  (north  of 
Trenton)  the  drumlins  are  very  large  and  some  outstanding  specimens 
exceed  200  feet  in  height  with  25  per  cent  slopes.    The  whole  area  is 
accompanied  by  wet  swampy  hollows.    Land  use  is  mainly  for  animal 
grazing  and  hay  and  oats  are  grown. 

The  land  use  and  other  information  are  presented  in  Table  1 
and  also  shown  in  Figures  2  and  3.     Industrial  and  urban  development 
is  small,  and  as  this  analysis  shows,  represents  an  insignificant 
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proportion  of  total  discharge  of  phosphorus  from  the  watersheds  to 
Lake  Ontario.    Forest  and  pasture  land  are  by  far  the  most  pre- 
dominant land  uses  in  the  area  of  the  four  watersheds.    The  largest 
urban  area  is  Belleville  with  6.4  square  miles  in  the  Moira  River 
basin  „    Napanee  is  the  only  municipality  served  by  a  water  pollution 
control  plant.    The  discharge  of  the  treatment  plant  enters  the  Napanee 
River  just  prior  to  confluence  with  Lake  Ontario. 

METHOD  OF  ANALYSIS  OF  PHOSPHORUS  DATA 
The  water  quality  data  were  collected  near  the  mouth  of  the 
four  streams  on  a  monthly  or  bimonthly  basis  during  the  study  period 
(October  1964  to  September  1967).    Streamflow  Information  was  ob- 
tained for  the  days  of  sampling  from  the  Canada  Department  of  Energy, 
Mines  and  Resources.    On  an  average  about  twenty- four  samples  were 
collected  on  each  stream  during  the  study  period.    The  samples  were 
obtained  on  random  days,  and  thus  represent  a  random  series. 

For  all  the  samples  collected  the  total  phosphorus  concentration 
level  (mg/1),  was  measured,  using  standard  methods  (A.P.H.A.  1960) 
(1).    In  Table  2,  the  annual  yields  based  on  mean  concentrations  and 
mean  flows  for  the  samples  collected  during  each  water  year  (October  1 
to  September  30)  are  presented. 


5. 

Figures  4,  5,  6  and  7  show  the  variations  of  the  phosphorus 
concentration  levels  on  sampling  days  as  well  as  the  flows  (ft  /m^-day) 
and  yields  of  phosphorus  (lbs/rrr-day)  for  the  four  streams .    These 
graphs  give     an  indication  of  the  relationship  between  the  three 
variables,  if  present.     In  any  attempt  to  correlate  these  variables 
statistically,  however,  one  must  keep  in  mind  that  the  yield  is  not 
an  independent  variable  since  it  is  a  product  of  the  concentration  and 
flow. 

In  Figure  8  the  average  monthly  yields  were  plotted  for  each 
of  the  four  streams.     Mean  monthly  yields  (lbs/month)  are  based  on 
mean  monthly  concentrations  and  the  mean  monthly  flows  (ftvday)  for 
the  study  period.     Before  the  regression  analysis  was  tried  for  different 
variables,   simple  log  relationships  between  concentration  and  flows 
were  plotted,  however,  these  graphs  are  not  presented  here  as  they  did 
not  show  any  significant  relationship. 

In  Figure  9,  the  ranges  between  yearly  maximum  and  minimum 
concentration  and  mean  values  have  been  plotted  for  the  four  streams. 
Figure  10  shows  the  range  between  maximum  and  minimum  yields. 


DISCUSSION 
A  number  of  researchers  in  the  past  have  drawn  attention  to  the 
importance  of  land  drainage  as  a  significant  source  of  phosphorus  to 


surface  waters.    In  1966,  Owen  and  Johnson  (2)  showed  a  statistically 
significant  relationship  between  the  yield  of  phosphorus  (lbs/mi2-day) 
with  flows  (ftVmi^-day)  and  developed  an  equation  which  allowed  fore- 
casting of  phosphorus  yield  for  streams .    However,  the  relationship  is 
somewhat  questionable  since  yield  and  flow  are  not  independent  vari- 
ables.   For  a  regression  analysis,  it  is  essential  that  random  samples 
be  obtained  from  the  values  of  the  dependent  variable  (y)  for  those  values 
of  the  independent  variable  (x)  included  in  the  study.    Also  when  (y)  is 
a  multiple  of  x  then  variation  in  x  will  not  explain  the  true  variation  in 
y  because  the  effect  of  the  variable  x  has  been  included  in  the  variable 
y.    The  yield  which  is  a  product  of  concentration  and  flow  cannot  be 
correlated  with  flows. 

Figures  4  to  7  show  that  total  phosphorus  content  (mg/1) 
fluctuated  markedly  in  some  months.    In  general,  the  fluctuations  did 
not  associate  with  changes  in  flow  (ft  /mi^-day)  and  appeared  to 
indicate  no  direct  or  inverse  relationship  between  phosphorus  concen- 
trations and  flows.    This  relationship  was  verified  by  calculating  the 
coefficient  of  correlation  between  the  two  variables  and  the  required 
critical  values  for  the  coefficient  of  correlation  at  the  5  per  cent  level 
of  significance  (Table  3).    By  comparing  the  calculated  and  required 
values  in  this  table  it  is  apparent  that  at  that  significance  level,  no 


relationship  exists  between  concentration  levels  and  flows.    This 
appears  to  indicate  that  streamflow  is  not  an  important  factor  for 
forecasting  the  phosphorus  concentrations  in  these  four  streams. 

An  attempt  was  also  made  to  correlate  the  concentration 
with  temperature  (°C)  of  the  sample  water,  and  the  dissolved  oxygen 
(DO)  (mg/1)  in  the  sample.    This  was  done  by  using  the  technique 
of  stepwise  multiple  regression  analysis.    However,  no  significant 
statistical  relationship  was  found,  and  the  equation  is  therefore  not 
included  in  this  report. 

Figure  8  shows  the  seasonal  variation  of  the  mean  monthly 
yield  (PO4  lbs/mi  -day)  of  phosphorus  for  each  watershed.    All  four 
watersheds  show  great  fluctuations  in  phosphorus  yields  between 
periods  of  high  and  low  streamflows.    About  18  to  25%  of  the  annual 
yield  was  discharged  from  these  streams  during  the  two  months  of 
December  and  January.    About  40  to  50%  of  the  annual  yield  was  dis- 
charged during  February,  March  and  April  as  compared  with  only  3  to 
22%  for  the  period  of  May  to  September.    By  comparison,  Engelbrecht 
and  Morgan  (3)  found  that  only  approximately  6%  of  the  total  annual 
yield  from  a  predominantly  agriculture  watershed  was  discharged 
during  a  period  of  April  to  September.     Some  seasonal  variation  in 
yield  is  intermediate  while  some  is  high. 
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Sawyer  (4),  in  his  study  of  Madison  Lake,  describes  a  mean 
annual  yield  of  phosphorus  from  agriculture  drainage  of  this  area 
equivalent  to  about  700  P04  lbs/mi2-yr.     Smith  (5)  found  a  yield  of 
114  P04  lbs/mi2  from  agriculture  drainage  during  a  period  of  April 
through  December.    Engelbrecht  and  Morgan  (6)  estimated  that  ap- 
proximately 665  P04  lbs/mi2-yr.  was  of  land  drainage  origin  on  an 
Illinois  water-shed  of  4,220  square  miles.     In  the  study  of  Eastern 
Lake  Ontario  Region,  the  mean  yield  of  phosphorus  as  P04  was  found 
to  be  332  lbs/mi2-yr.  from  the  Moira,   237  lbs/mi2-yr.  from  the  Salmon, 
810  lbs/mi2-yr.  from  the  Napanee  and  375  lbs/mi  -yr.  from  Wilton 

Creek. 

The  four  watersheds  studied  here  have  little  industry  and 
urban  area.    The  Napanee  River  receives  the  effluent  from  the  sewage 
treatment  plant  serving  the  Town  of  Napanee  which  increased  the  yield 
of  phosphorus  for  this  watershed  to  810  lbs  per  square  mile  per  year  (lb/ 
mi2-yr).  With  the  phosphate  input  from  the  Napanee  water  pollution 
control  plant  deducted,  the  phosphorus  yield  from  this  watershed 
reduces  to  654  lbs/mi2- yr.    No  attempt  was  made  to  investigate  in 
detail  the  reason  for  having  high  yield  at  Napanee,  however,  a  visual 
observation  and  discussion  with  local  people  Indicates  that  this  high 
value  may  be  due  to  the  flooding  of  lands  and  existence  of  Depot  lakes 
on  the  river,  and  also  because  of  effluent  coming  from  paper  industry 
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into  the  stream.    From  the  knowledge  that  land  use  in  the  four  streams 
is  predominantly  forest  and  grazing ,  it  can  be  apparently  concluded 
that  drainage  from  these  two  land  uses  provides  the  main  source  of 
phosphorus  to  the  rivers  in  this  region  throughout  the  year. 

The  bar  diagrams  in  Figure  9  indicate  the  maximum,  minimum , 
average  concentration  values  based  on  study  period  (1964-67) .    While 
the  available  data  does  not  allow  to  comment  on  developing  trends  at 
this  time,  it  Is  evident  that  with  more  data  such  bar  diagrams  are  use- 
ful in  establishing  long  term  trends  in  the  variations  of  the  yield  of 
phosphorus  from  the  area  under  study  or  elsewhere. 

The  equation  obtained  from  the  multiple  regression  analysis 
relating  phosphorus  yield  to  major  land  uses  on  four  streams  is  shown 
in  Table  3.    The  land  uses  for  this  region  are  forest,  horticulture, 
urban,   pasture  and  grazing  land.    The  cursory  examination  appears 
to  indicate  a  relationship  between  phosphorus  concentration  and  horti- 
cultural and  urban  land  uses  on  the  four  watersheds.    However,  as 
only  four  streams  were  considered,  the  number  of  observations  for 
this  multiple  regression  was  limited  to  four,  and  any  attempt  at 
regression  analysis  having  more  variables  than  observations  would 
be  meaningless,  (the  number  of  degrees  of  freedom  would  be  negative). 
Therefore,  after  the  addition  of  two  variables  the  multiple  regression 
was  terminated.    Although  the  multiple  correlation  coefficient  is  high, 
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r  =  0.959,  It  is  below  the  critical  r    =  0.999  required  at  a  5  per  cent 
level  of  significance.    It  is,  therefore,  concluded  that  with  the  data 
used  no  significant  relationship  can  be  established  between  phos- 
phorus yield  and  land  uses  and  the  established  forecasting  equation 
does  not  produce  reliable  results.    The  same  holds  true  for  other 
equations  based  on  simple  regression  analysis  for  the  four  streams 
(Table  3). 

SUMMARY  AND  CONCLUSION 

Concentration  levels  of  total  phosphorus  in  four  major  streams 
of  eastern  Lake  Ontario  region  were  investigated.    The  watershed  of 
all  these  streams  is  predominantly  wooded  and  grazing  land.    There 
is  neither  much  organized  area  nor  much  industry  in  these  watersheds, 
with  the  result  that  the  contribution  of  phosphorus  would  originate 
mainly  from  forest  and  pasture  lands.    The  finding  of  this  report  shows 
a  significant  yield  of  phosphorus  from  the  non-urbanized  areas. 

No  relationship  was  found  to  exist  between  level  of  concen- 
tration in  the  tributary  drainage  and  the  rate  of  streamflow.    The  yield, 
however,  depended  on  the  fluctuations  of  flows  because  total  yield 
incorporates  within  itself  the  effect  of  flows.    Pronounced  seasonal 
variations  were  noted  in  the  yields  of  phosphorus  from  forest  and 
grazing  land  areas.    About  50%  of  the  annual  yield  occurred  in  February, 
March  and  April . 
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It  can  be  concluded  from  this  study  that  there  appears  to  be 
no  relationship  between  the  level  of  concentration  of  phosphorus  and 
streamflows.    However,  it  was  found  that  concentrations  of  phosphorus 
in  four  streams  was  about  twice  as  high  in  August  than  in  March „    This 
appears  to  differ  from  values  reported  by  Owen  and  Johnson  (2). 

The  analysis  presented  here  could  in  general  be  used  for  the 
assessment  of  yields  from  other  streams  where  adequate  water  quality 
monitoring  data  are  available.    But  in  order  to  obtain  statistically 
significant  results,  water  quality  data  for  longer  periods  should  be 
available.    To  increase  the  reliability  of  the  findings  presented  in 
this  case  study,  it  is  desirable  to  have  more  data,' about  15  years. 
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TABLE  1  -    CHARACTERISTICS  OF  FOUR  SELECTED  WATERSHEDS 
IN  EASTERN  LAKE  ONTARIO  REGION 


Watersheds 


Moira 
River 


Salmon 
River 


Napanee 
River 


Wilton 
Creek 


Area  Estimated  Pop. 

(sq.ml.)     Rural         Urban 


Major 

Land  Use 

(sq.ml,) 


Total  Effluent  Dis- 
charge above  sample 
point  (IMGD) 


1,040 


344 


300 


10,000     3,600 


10,000 


11,00        7,600 


4.34         3,000 


urban  9 . 1 
woodland  693 
pasture  & 
forage 
land  261.6 

urban  3 .68 
woodland  203 
pasture  & 
forage 
land  112.34 

urban  3 
woodland  165 
pasture  & 
forage 
land  129.7 

urban  .5 
woodland  10.' 
pasture  & 
forage 
land  38.1 


M  -      -6 

I    -     2.7 


I   -       .025 
(seasonal 
operation) 


M  -      .9 
I    -      .9 


I    -      0.03 


No.  of         Municipal 

Indus-         Treatment 

tries  (1)      Plants  (M) 


Mean  Discharge 
(1000  ft3/day) 
During  Study  Pd. 
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85102.6 


28624.9 


24935.0 


3647.4 


oo 


TABLE  2  -  MEAN  CONCENTRATION  LEVEL,   MEAN  FLOW  AND  ANNUAL  YIELD  OF  PHOSPHORUS 
FROM  FOUR  WATERSHEDS  IN  THE  EASTERN  LAKE  ONTARIO  REGION 


Moira  River 

Mean  Concentration 
(mq/1) 

Mean  Flow 
(cfs/mi^) 

Avg.  Annual  Yiel 
(lbs/mi2-  yr) 

d 

1964-65 

1965-66 

1966-67 

1964-65 

1965-66 

1966-67 

1964-65     1965-66 

1966-67 

.22 

.19 

.138 

.7165 

1.048 

1.077 

310.7         392.4 

292.9 

Salmon  River 

.22 

.09 

.087 

.7994 

1.052 

1.038 

346.5         186.6 

177.9 

Napanee  River 

.74 

.42 

.195 

.786 

1.073 

1.027 

1146.3         888.1 

394.6 

Wilton  Creek 

.24 

.17 

.189 

.7823e 

1.099 

1.037 

369.7         368.1 

386.1 

TABLE  3  -  COEFFICIENT  OF  CORRELATION  BETWEEN  PHOSPHORUS  CONCENTRATION  (mg/l),y, 
AND  FLOW  (ft3/m2/day) ,  x,  FOR  FOUR  WATERSHEDS  IN  THE  EASTERN  LAKE  ONTARIO  REGION 


■ 

Watershed 

Estimating  Equation 

Computed  Correlation 
Coefficient  (r) 

Critl 
Level 

ca 
of 

1  Value  of  r  at  5% 
Significance  Req'd 

Moira  River 

y  =  .191  -   c024x 

.292 

.423 

Salmon  River 

y  =  .146  -   .028x 

.262 

.367 

Napanee  River 

y  =  .459  -   .121x 

.439 

.367 

Wilton  Creek 

y  -  .216  -   .041x 

.269 

.396 

Whole  Region 

y  =  -.78  +  .832x1  +  . 

296x 

i 

.959 

.999 

where 

y      - 

X2" 

is 
is 
is 

Phosphorus  Concentration 
urban ,  and 
horticulture  areas 

tn 


POPULATION    CENTRES 
AREA   WHERE  LAND  USE    IS 
3  PREDOMINANTLY  WOODLAND 
AREA  WHERE    LAND   USE    IS 
CROPLAND   AND   PASTURE 


LAKE    ONTARIO 


FIG.  2 
LAND  USE 
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FIG.     4 

MOIRA      RIVER 

VARIATION      OF      PHOS,    FLOW,      YIELD 

DURING      STUDY      PERIOD      1964  -  1967 
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FIG. 

SALMON       RIVER 
VARIATION      OF      PHOS,     FLOW,      YIELD 
DURING      STUDY       PERIOD 
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FIG      6 

NAPANEE      RIVER 
VARIATION      OF      PMOS,     FLOW,     YIELD 
DURING      STUDY      PERIOD     1964  - 1967 
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FIG.     7 

WILTON       CREEK 
VARIATION      OF      PHOS,     FLOW,     YIELD 
DURING      STUDY      PERIOD      1964  -  1967 
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FIG      t 

COMPARISON       OF       MEAN       MONTHLY 
PHOSPHOROUS       YIELD       FOR 
***\         SELECTED      WATERSHEDS 
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FIG.  9 
YEARLY     MAXIMUM  ,  MINIMUM  ,  AVERAGE     CONCENTRATION    OF    TOTAL     PHOSPHOROUS      (964  -  67 
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